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ABSTRACT 



J> i We present an improved comparison of the strong coupling of the gluon to light {qi = u+d+s), 

CN ' c, and b quarks, determined from multijet rates in flavor-tagged samples of hadronic 

o : 

, decays recorded with the SLC Large Detector at the SLAC Linear Collider between 1993 

00 ' and 1995. Flavor separation among primary qiTji, cc, and bb final states was made on the 



basis of the reconstructed mass of long-lived heavy-hadron decay vertices, yielding tags 
with high purity and low bias against > 3-jet final states. We find: a'^/ag'^^ = 1.036 it 
O : 0.043 {stat.)tl-'oil isyst.)tlfil (theory) and a'jaf' = 1.004±0.018 {staQt'ofsl {sysQt'of^l 



(theory). 

Submitted to Physical Review D 
*Work supported by Department of Energy contracts: DE-FG02-91ER40676 (BU), DE-FG03-91ER40618 
(UCSB), DE-FG03-92ER40689 (UCSC), DE-FG03-93ER40788 (CSU), DE-FG02-91ER40672 (Col- 
orado), DE-FG02-91ER40677 (Illinois), DE-AC03-76SF00098 (LBL), DE-FG02-92ER40715 (Mas- 
sachusetts), DE-FC02-94ER40818 (MIT), DE-FG03-96ER40969 (Oregon), DE-AC03-76SF00515 (SLAC), 
DE-FG05-91ER40627 (Tennessee), DE-FG02-95ER40896 (Wisconsin), DE-FG02-92ER40704 (Yale); 
National Science Foundation grants: PHY-91-13428 (UCSC), PHY-89-21320 (Columbia), PHY-92- 
04239 (Cincinnati), PHY-95-10439 (Rutgers), PHY-88-19316 (Vanderbilt), PHY-92-03212 (Washing- 
ton); The UK Particle Physics and Astronomy Research Council (Brunei, Oxford and RAL); The 
Istituto Nazionale di Fisica Nucleare of Italy (Bologna, Ferrara, Frascati, Pisa, Padova, Perugia); 
The Japan-US Cooperative Research Project on High Energy Physics (Nagoya, Tohoku); The Korea 
Research Foundation (Soongsil, 1997). 



1. Introduction 



In order for Quantum Chromodynamics (QCD) |]T| to be a gauge-invariant renormalis- 
able field theory it is required that the strong couphng between quarks (g) and gluons 
((?), cts, be independent of quark flavor. This basic ansatz can be tested directly in e+e~ 
annihilation by measuring the strong coupling in events of the type e+e~ — > qqg for 
specific quark flavors. Whereas an absolute determination of as using such a technique 
is limited, primarily by large theoretical uncertainties, to the 5%-level of precision p[, 
a much more precise test of the flavor-independence can be made from the ratio of the 
couplings for different quark flavors, in which most experimental errors and theoreti- 
cal uncertainties cancel. Furthermore, the emission of gluon radiation in bb events is 
expected to be modified relative to that in qiql {qi =u+d+s) events due to the large 
6-quark mass, and comparison of the rates for —>■ bbg and QiWiQ "^^Y allow 

measurement of the running mas^ of the 6-quark, mb(M^o)[|. Finally, in addition to 
providing a powerful test of QCD, such measurements allow constraints to be placed on 
physics beyond the Standard Model. For example, a flavor- dependent anomalous quark 
chromomagnetic moment would modify the emission rate of gluons for the different 
quark flavors, and would manifest itself in the form of an apparently flavor-dependent 
strong coupling. 

The first such comparisons, of as for c or 6 quarks with as for all flavors, were made 

at the PETRA e"'"e~ collider at cm. energies in the range 35 < ^/s < 47 GeV and were 

limited in precision to da^/af^ = 0.41 and da^af^ = 0.57 [0 due to the small data 

sample and hmited heavy-quark tagging capability. These studies made the simplifying 

assumptions that a^ = a^'^^ and a^ = a"'^'', respectively. More recently, measurements 
*Use of the modified minimal subtraction renormalisation scheme is implied throughout this 
paper. 

^The DELPHI Collaboration has recently measured the three-jet rate ratio R^/K^'^^ to a precision 
of ±0.009, and, under the assumption of a flavor-independent strong coupling, derived a value of the 
running 6-mass this issue will be discussed in Section 6. 



2 



made at the resonance have benefitted from the use of micro-vertex detectors 
for improved heavy-quark tagging. Samples of tagged bb events recorded at LEP have 
been used to test flavor-independence to a precision of 6al/af = 0.012 [g i, but these 
measurements were insensitive to any differences among values for the non-6-quarks. 
The ALEPH Collaboration also measured a'l'^/a^'^^ to a precision of ±0.023 but in 
this case there is no sensitivity to a different as for c and b quarks. 

The OPAL Collaboration has measured a{/af^ for all five flavors / with no as- 
sumption on the relative value of for the different flavors |]TD[, and has verified 
fiavor-independence to a precision of dal/af- = 0.026, da'^/af' = 0.09, 5al/af^ = 
0.15, = 0.20, and 6a^/af'' = 0.21. In that analysis the precision of the test 

was limited by the kinematic signatures used to tag c and light-quark events, which 
suffer from low efficiency and strong biases against events containing hard gluon radi- 
ation. In our previous study |[Tl|] we used hadron lifetime information as a basis for 
separation of bb, cc and light-quark events with relatively small bias against 3-jet final 
states. We verified fiavor-independence to a precision of Sa^/af'' = 0.06, 6a'^/af'' = 
0.17, and Sa^^Jaf = 0.04. 

Here we present an improved test of the fiavor-independence of strong interactions 
using a sample of hadronic decay events produced by the SLAC Linear Collider 
(SLC) and recorded in the SLC Large Detector (SLD) in data-taking runs between 1993 
and 1995. The precise tracking capability of the Central Drift Chamber and the 120- 
million-pixel CCD-based Vertex Detector (VXD2), combined with the stable, micron- 
sized beam interaction point (IP), allowed us to reconstruct topologically secondary 
vertices from heavy-hadron decays with high efficiency. High-purity samples of Z^ 

bb{g) and Z° — > cc(g) events were then tagged on the basis of the reconstructed mass 
and momentum of the secondary vertex. Events containing no secondary vertex and 
no tracks significantly displaced from the IP were tagged as a high-purity Z^ QiQiio) 
event sample. The method makes no assumptions about the relative values of a^, 
and a"'^*. Furthermore, an important advantage of the method is that it has low bias 
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against > 3-jet events. In addition to using an improved flavor-tagging technique, 
this analysis utilises a data sample three times larger than that used for our previous 
measurement, and allows us to test the flavor independence of strong interactions to a 
precision higher by roughly a factor of three. Finally, quark mass effects in —* qqg 
events have recently been calculated |l3l at next-to-leading order in perturbative 
QCD, and are non-negligible on the scale of our experimental errors; we have utilised 
these calculations in this analysis. 



2. Apparatus and Hadronic Event Selection 

This analysis is based on roughly 150,000 hadronic events produced in e"'"e~ annihilations 
at a mean center-of-mass energy of ^/s = 91.28 GeV. A general description of the SLD 
can be found elsewhere [|1^]. The trigger and initial selection criteria for hadronic 
decays are described in Ref. ||T3|. This analysis used charged tracks measured in the 
Central Drift Chamber (CDC) pl and in the Vertex Detector (VXD2) 0. Momen- 
tum measurement is provided by a uniform axial magnetic field of 0.6T. The CDC 
and VXD2 give a momentum resolution of (Jp^/pi_ = 0.01 © 0.0026p_L, where p± is the 
track momentum transverse to the beam axis in GeV/c. In the plane normal to the 
beamline the centroid of the micron-sized SLC IP was reconstructed from tracks in 
sets of approximately thirty sequential hadronic Z^ decays to a precision of ct/p ^ 7 
/im. Including the uncertainty on the IP position, the resolution on the charged-track 
impact parameter [d) projected in the plane perpendicular to the beamline is = 
11©70/ {p± sin^/^ 6) /im, where 6 is the track polar angle with respect to the beamline. 



The event thrust axis |]T8[ was calculated using energy clusters measured in the Liquid 
Argon Calorimeter [|19| . 

A set of cuts was applied to the data to select well-measured tracks and events 
well contained within the detector acceptance. Charged tracks were required to have a 
distance of closest approach transverse to the beam axis within 5 cm, and within 10 cm 
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along the axis from the measured IP, as well as | cos^| < 0.80, and p± > 0.15 GeV/c. 
Events were required to have a minimum of seven such tracks, a thrust axis polar 
angle w.r.t. the beamline, 6t, within |cos^t| < 0.71, and a charged visible energy 
Eyis of at least 20 GeV, which was calculated from the selected tracks assigned the 
charged pion mass. The efficiency for selecting a well-contained Z° qqig) event was 
estimated to be above 96% independent of quark flavor. The selected sample comprised 
77,896 events, with an estimated 0.10 ± 0.05% background contribution dominated by 
T^T^ events. 

For the purpose of estimating the efficiency and purity of the event flavor-tagging 
procedure we made use of a detailed Monte Carlo (MC) simulation of the detector. The 
JETSET 7.4 PD[ event generator was used, with parameter values tuned to hadronic 



e"'"e annihilation data , combined with a simulation of 5-hadron decays tuned p2 



to T(4S') data and a simulation of the SLD based on GEANT 3.21 ||23|. Inclusive 
distributions of single-particle and event-topology observables in hadronic events were 
found to be well described by the simulation ||l5l . Uncertainties in the simulation were 
taken into account in the systematic errors (Section 5). 



3. Flavor Tagging 

Separation of the accepted event sample into tagged flavor subsamples was based on the 
invariant mass of topologically-reconstructed long-lived heavy-hadron decay vertices, 
as well as on charged-track impact parameters in the plane normal to the beamline. 
In each event a jet structure was deflned as a basis for flavor-tagging by applying the 
'JADE' jet-flnding algorithm to the selected tracks; a value of the normalised jet-jet 
invariant-mass parameter Uc = 0.02 was used. The impact parameter of each track, d, 
was given a positive (negative) sign according to whether the point-of-closest approach 
to its jet axis was on the same side (opposite side) of the IP as the jet. Charged tracks 
used for the subsequent event flavor-tagging were further required to have at least 40 
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hits in the CDC, with the first hit at a radial distance of less than 39 cm from the 
beamhne, at least one VXD2 hit, a combined CDC + VXD2 track fit quality of Xdof 
< 5, momentum p > 0.5 GeV/c, a distance of closest approach transverse to the beam 
axis within 0.3 cm, and within 1.5 cm along the axis from the measured IP, and an 
error on the impact parameter, cr^, less than 250/im. Tracks from identified K° and A 
decays and 7 conversions were removed. 

In each jet we then searched for a secondary vertex (SV), namely a vertex spatially 
separated from the measured IP. In the search those tracks were considered that were 
assigned to the jet by the jet-finder. Individual track probability-density functions in 
3-dimensional co-ordinate space were examined and a candidate SV was defined by 
a region of high track overlap density; the method is described in detail in [^]. A 
SV was required to contain two or more tracks, and to be separated from the IP by 
at least 1 mm. We found 14,096 events containing a SV in only one jet, 5817 events 
containing a SV in two jets, and 54 events containing a SV in more than two jets. The 
selected SVs comprise, on average, 3.0 tracks. These requirements preferentially select 
SVs that originate from the decay of particles with relatively long lifetime. In our 
simulated event sample a SV was found in 50% of all true 6-quark hemispheres, in 15% 
of true c-quark, and in < 1% of true light-quark hemispheres p5[, where hemispheres 
were defined by the plane normal to the thrust axis that contains the IP. 

Due to the cascade structure of B-hadron decays, not all the tracks in the decay 
chain will necessarily originate from a common decay point, and in such cases the SV 
may not be fully reconstructed in bb events. Therefore, we improved our estimate of 
the SV by allowing the possibility of attaching additional tracks. First, we defined 
the vertex axis to be the straight line joining the IP and the SV centroids, and D 
to be the distance along this axis between the IP and the SV. For each track in the 
jet not included in the SV the point of closest approach (POCA), and corresponding 
distance of closest approach, T, to the vertex axis were determined. The length, L, of 
the projection of the vector joining the IP and the POCA, along the vertex axis was 
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then calculated. Tracks with T < 1.0 mm, L > 0.8 mm and L/D > 0.22 were then 
attached to the SV. On average 0.5 tracks per SV were attached in this fashion. 

The invariant mass, Mch, of each SV was then calculated by assigning each track 
the charged pion mass. In order to account partially for the effect of neutral particles 
missing from the SV we apphed a kinematic correction to the calculated Mcu- We added 
the momentum vectors of all tracks forming the SV to obtain the vertex momentum, 

— * 

Pvtxi and evaluated the magnitude of the component of the vertex momentum tranverse 
to the vertex axis, Pj. In order to reduce the effect of the IP and SV measurement 
errors, the vertex axis was varied within an envelope defined by all possible cotangents 
to the error eUipsoids of both the IP and the SV, and the minimum Pt was chosen. We 



The distributions of M^t^ and P^tx are shown in Fig. 1; the data are reproduced by 
the simulation, in which the primary event-flavor breakdown is indicated. The region 
Mytx > 2 GeV/c^ is populated predominantly by Z° bb events, whereas the region 
Mytx < 2 GeV/c^ is populated roughly equally by bb and non-66 events. 

In order to optimise the separation among flavors we examined the two-dimensional 
distribution of P^tx vs. M^tx- The distribution for events containing a SV is shown in 
Fig. 2 for the data and simulated samples; the data (Fig. 2a) are reproduced by the 
simulation (Fig. 2b). The distributions for the simulated subsamples corresponding to 
true primary bb, cc, and qiql events are shown in Figs. 2c, 2d and 2e respectively. 

In order to separate bb and cc events from each other, and from the qiqi events, 
we defined the regions: (A) M^tx > 1-8 © P^tx + 10 < IbM^tx] (B) M^^ < 1-8 © 
P^t^ > 5 © P^tx + 10 > IbM^tx] where M^,tx {Pvtx) is in units of GeV/c^ (GeV/c); 
(G) all remaining events containing a SV. The boundaries of regions (A) and (B) are 
indicated in Figs. 2c and 2d, respectively, and all three regions are labelled in Fig. 2f. 
The 6-tagged sample (subsample 1) was defined to comprise those events containing 
any vertex in region (A). For the remaining events containing any vertex in region 
(B) we examined the distribution of the impact parameter of the vector P^tx w.r.t. 



then defined the Pj-corrected vertex mass, M^t^ — \jM\^ -|- P^ -|- \Pt 
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the IP, 5ytx (Fig. 3); according to the simulation true primary cc events dominate the 
population in the region d^tx < 0.02 cm. Therefore, we defined the c-tagged sample 
(subsample 2) to comprise those events in region (B) with S^tx < 0.02 cm. 

Events containing no selected SV were then examined. For such events the dis- 
tribution of Nsig, the number of tracks per event that miss the IP hj d > 2ad, is 
shown in Fig. 4. The ucis-tagged sample (subsample 3) was defined to comprise those 
events with Nsig = 0. All events not assigned to subsamples 1,2 or 3 were defined to 
comprise the untagged sample (subsample 4). Using the simulation we estimated that 
the efficiencies e^'' for selecting events (after acceptance cuts) of type i {i = b, c, uds, ) 
into subsample j (1 < j < 4), and the fractions of events of type i in subsample 
j, are {e,nY^ = (61.5 ± 0.1%, 95.5 ± 0.1%), {e, Uf" = (19.1 ± 0.1%, 64.4 ± 0.3%) 
and (e, n)^"'^'' = (56.4 ± 0.1%, 90.6 ± 0.1%). The composition of the untagged sam- 
ple (subsample 4) was estimated to be U^""^' = 59.3 ± 0.1%, U^'' = 24.1 ± 0.1% and 
= 16.6 ± 0.1%. The errors on these values are discussed in Section 5. 



4. Jet Finding 

For the study of fiavor-independence the jet structure of events was reconstructed in 
turn using six iterative clustering algorithms. We used the 'E', 'EO', 'P', and 'PO' 
variations of the JADE algorithm, as well as the 'Durham' ('D') and 'Geneva' ('G') 
algorithms ^6]. In each case events were divided into two categories: those containing 
(i) two jets, and (ii) three or more jets. The fraction of the event sample in category (ii) 
was defined as the 3-jet rate R^. This quantity is infrared- and coUinear-safe and has 



been calculated to 0{al) in perturbative QCD ||26|, For each algorithm we repeated 
the subsequent analysis successively across a range of values of the normalised jet-jet 
invariant-mass parameter i/c, 0.005 < Uc ^ 0.12. The ensemble of results from the 
different Uc values was used to cross-check the consistency of the method. In the final 
stage an 'optimal' yc value was chosen for each algorithm so as to minimise the overall 
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error on the analysis, and the spread in results over the algorithms was used to assign 
an additional uncertainty (Section 7). 

Each of the six jet-finding algorithms was applied to each tagged-event subsample j, 
1 < J < 3 (Section 3), as well as to the global sample of all accepted events {'alP). For 
each algorithm the 3- jet rate in each subsample was calculated, and the ratios Rl/Rf'', 
in which many systematic errors should cancel, were then derived. As an example the 
Ry Rf" are shown as a function of yc for the JADE EO algorithm in Fig. 5a. The 
results of the corresponding analysis applied to the simulated event sample are also 
shown; the simulation reproduces the data. Similar results were obtained for the other 
jet algorithms (not shown). 

For each algorithm and i/c value the R^ for each of the i quark types {i — b, c, uds) 
was extracted from a simultaneous maximum likelihood fit to 77-2 and 77-3, the number of 
2-jet and 3-jet events, respectively, in the fiavor-tagged subsample (1 < j < 3), using 
the relations: 

i=uds,c,b 

4 ^ E (4-3)^3 + ^;l3)(i - ^3)) fN . (1) 

i=uds,c,b 

Here is the total number of events after correction for the event selection efficiency, 
and /* is the Standard Model fractional hadronic width for decays to quark type i. 
The T/c-dependent 3x3 matrices £(2-^^2) ^(3-^3) efficiencies for an event of 

type i, with 2- or 3-jets at the parton level, to pass all cuts and enter subsample j as 
a 2- or 3-jet event, respectively. Similarly, the 3x3 matrices £(2^3) and £(3^2) 
the efficiencies for an event of type i, with 2- or 3-jets at the parton level, to pass all 
cuts and enter subsample j as a 3- or 2-jet event, respectively These matrices were 
calculated from the Monte Carlo simulation, and the systematic errors on the values 
of the matrix elements are discussed in Sections 5 and 6. 

This formalism explicitly accounts for modifications of the parton-level 3-jet rate 
due to hadronisation, detector effects, and flavor-tagging bias. The latter effect is 

9 



evident, for the EO algorithm, in Fig. 5a, where it can be seen that the measured 
values of Rl/Rf^ are below unity for subsamples j = 1,2 and 3, implying that the 
flavor tags preferentially select 2-jet rather than 3-jet events. For example, at yc = 
0.02 the normalised difference in efficiencies for correctly tagging a 2-jet event and 
a 3-jet event of type i in subsample j are i3^'^=5.7%, i3^'''=14.5%, and B^''"^^=4.1%, 
where B^"" = {e2'^2 ~ ^i^-^s) / ^2-^2'^ these biases are considerably smaller than those found 



in [|T^, which resulted from the kinematic signatures employed for flavor-tagging. It 
should be noted that, as a corollary, the untagged event sample, subsample 4, contains 
an excess of 3-jet events (Fig. 5a). Similar results were obtained for the other jet 
algorithms (not shown). 

Equations 1 were solved using 2- and 3-jet events defined in turn by each of the 
six jet algorithms to obtain the true 3-jet rates in qiql, cc and hb events, Kif''^, 

-R3 and R\ respectively. Redefining Rf^ = T,b^c,udsf^^h unfolded ratios R^'^^ /Rf\ 
R^Rf^ and R\/ Rf^ are shown in Fig. 5b for comparison with the raw measured values 
shown in Fig. 5a. 

For the test of the flavor-independence of strong interactions it is more convenient to 



uds 
3 



consider the ratios of the 3-jet rates in heavy- and light-quark events, namely R^/R; 
and R'^/R^'^^. These were derived from the unfolded R^'^^, R^ and -R3 values, and the 
systematic errors on the ratios are considered in the next sections. 



5. Experimental Systematic Errors 

We considered sources of experimental systematic uncertainty that potentially affect 
our measurements of R"^/ R^'^'^ and R\/ R^'^^ . These may be divided into uncertain- 
ties in modelling the detector and uncertainties on experimental measurements serving 
as input parameters to the underlying physics modelling. In each case the error was 
evaluated by varying the appropriate parameter in the Monte Carlo simulation, recal- 
culating the matrices e, performing a new fit of Eq. 1 to the data, rederiving values 
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of i^g/i^g'^* and R\/K^'^^, and taking the respective difference in results relative to our 
standard procedure as the systematic uncertainty. 

In the category of detector modelling uncertainty we considered the charged-particle 
tracking efficiency of the detector, as well as the smearing applied to the simulated 
charged-particle impact parameters in order to make the distributions agree with the 
data. An extra tracking inefficiency of roughly 3.5% was applied in the simulation in 
order to make the average number of charged tracks used for flavor-tagging agree with 
the data. We repeated the analysis in turn without this efficiency correction, and with 
no impact-parameter smearing, in the simulation. 

A large number of measured quantities relating to the production and decay of 
charm and bottom hadrons are used as input to our simulation. In hh events we 
have considered the uncertainties on: the average charged multiplicity of 5-hadron 
decays, the i?-hadron fragmentation function, the production rate of 6-baryons, the B- 
meson and 5-baryon lifetimes, the inclusive production rate of mesons in i?-hadron 
decays, and the branching fraction for hh, f^. In cc events we have considered 

the uncertainties on: the branching fraction for cc, the charmed hadron 

fragmentation function, the inclusive production rate of mesons, and the charged 
multiplicity of charmed hadron decays. We also considered the rate of production of 
secondary hb and cc from gluon splitting in qqg events. The values of these quantities 
used in our simulation and the respective variations that we considered are listed in 
Table 1. 

Statistical errors resulting from the finite size of the Monte Carlo event sample were 
estimated by generating 1,000 toy Monte Carlo datasets of the same size as that used in 
our data correction procedure, evaluating the matrices e (Eq. 1) for each, unfolding the 
data, and calculating the r.m.s. deviation of the distributions of the resulting R^/ R^'^^ 
and R\/Rf' values. 

As an example, for the EO algorithm at Uc = 0.02 the errors on i^g/ K^'^'^ and R^/ R^'^'^ 
from the above sources are listed in Table [1|. The dominant physics contributions to 
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5R\/R^°-^ result from limited knowledge of the average S-hadron decay multiplicity 
and the S-hadron fragmentation function. The uncertainties in and in the charmed 
hadron fragmentation function produce the dominant variations in R^/ R^'^'^ . Con- 
tributions from i?-hadron lifetimes, the fraction of in B meson decays, &-baryon 
production rates, and the charm hadron decay multiplicity are relatively small. 

For each jet algorithm and yc value all of the errors were added in quadrature to 
obtain a total experimental systematic error on R\/R^'^^ and R\/ K^'^^. The choice of 
an optimal yc value is discussed in Section 6, and the combination of results from the 
six jet algorithms is discussed in Section 7. 



6. Theoretical Uncertainties and Translation to as Ratios 

We considered sources of theoretical uncertainty that potentially affect our measure- 
ments. The ratios R'^/R^'^^ and R\/R^'^^ derived in Section 4 were implicitly corrected 
for the effects of hadronisation, and we have estimated the uncertainty in this correc- 
tion. Furthermore, the> 3-jet rate in heavy-quark events is modified relative to that 
in light-quark events by the effect of the non-zero quark mass. This effect needs to be 
taken into account in the translation between the jet-rate ratios and the corresponding 
ratios of strong couplings ctg/a"'^^ and a^/a"'^*. We have used 0(0;^) calculations to 
perform the mass-dependent translation, and have estimated the related uncertain- 
ties due to the value of the 6-quark mass, as well as higher-order perturbative QCD 
contributions. 



6.1 Hadronisation Uncertainties 

The intrinsically non-perturbative process by which quarks and gluons fragment into 
the observed final-state hadrons cannot currently be calculated in QCD. Phenomeno- 
logical models of hadronisation have been developed over the past few decades and have 
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been implemented in Monte Carlo event-generator programs to facilitate comparison 
with experimental data. We have used the models implemented in the JETSET 7.4 
and HERWIG 5.9 p8| programs to study hadronisation effects; these models have been 



extensively studied and tuned to provide a good description of detailed properties of 
hadronic final states in e"'"e~ annihilation; for a review of studies at the resonance 
see [|^. Our standard simulation based on JETSET 7.4 was used to evaluate the 
efficiency and purity of the event-flavor tagging, as described in Section 4, as well as 
for the study of experimental systematic errors described in Section 5. 

We investigated hadronisation uncertainties by calculating from the Monte Carlo- 
generated event sample the ratios: 

_ I R\ \ , ( R\ \ 



J^uds j 'I J^uds j 

\ 3 / parton \ 3 / hadron 

where i = c or b, parton refers to the calculation of the quantity in brackets at the 
parton-level, and hadron refers to the corresponding hadron-level calculation using 
stable final-state particles. We recalculated these ratios by changing in turn the pa- 
rameters Qo and aq in the JETSET program^ and generating 1-million-event samples. 
We also recalculated these ratios by using the HERWIG 5.9 program with default pa- 
rameter settings. For each variation we evaluated the fractional deviation Ar^ w.r.t. 
the standard value: 

Ari 



fr' - Ti 



and the corresponding deviations on Ry R^'^'^. As an example, for the EO algorithm and 

He = 0.02 the deviations are hsted in Table 1. The deviations were added in quadrature 

to define the systematic error on R\/ R^'^^ due to hadronisation uncertainties. 

^Qo (GeV) controls the minimum virtual mass allowed for partons in the parton shower; we con- 
sidered a variation around the central value, 1.0, of ^o"- (GeV/c) is the width of the Gaussian 
distribution used to assign transverse momentum, w.r.t. the color field, to quarks and antiquarks 
produced in the fragmentation process; we considered a variation around the central value, 0.39, of 
±0.04. 
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6.2 Choice of y,. Values 



For each jet algorithm and yc value the statistical and experimental systematic errors 
and hadronisation uncertainty on each R\/K^'^'^ were added in quadrature. No strong 
dependence of this combined error on Uc was observed [0, but an 'optimal' Uc value for 
each algorithm was then identified that corresponded with the smallest error. In the 
case of the E and G algorithms slightly larger Uc values were chosen so as to ensure that 
the O(a^) calculations for massive quarks were reliable |31|. The chosen Uc value for 
each algorithm is listed in Table 2, together with the corresponding values of the ratios 
R'^/R^'^'^ and R'^/R^'^'^, as well as the statistical and experimental-systematic errors and 
hadronisation uncertainties. 



6.3 Translation to as Ratios 

The test of the flavor-independence of strong interactions can be expressed in terms 
of the ratios al/a^'^^ {i = c or b). Recalling that with our definition R3 is the rate 
of production of 3 or more jets, can be derived from the respective measured 

ratio R\/ R^'^^ using the next-to-leading-order perturbative QCD calculation: 

where the coefficients A, B and C represent, respectively, the leading-order (LO) per- 
turbative QCD coefficient for the 3-jet rate, the next-to-leading-order (NLO) coefficient 
for this rate, and the leading-order coefficient for the 4-jet rate. Next-to- leading-order 
contributions to the 4-jet rate, and contributions from > 5-jet rates, are represented by 
the terms of O(a^). These coefficients depend implicitly upon the jet algorithm as well 
as on the scaled-invariant-mass-squared jet resolution parameter for clarity these 
dependences have been omitted from the notation. For massless quarks calculations of 
the coefficients A, B and C have been available for many years ||2^, |27[| . 

For many observables at the pole the quark mass appears in terms proportional 
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to the ratio m'^/M^, and the effects of non-zero quark mass can be neglected. For 
the jet rates, however, mass effects can enter via terms proportional to m'^ / [ycM"^) . 
For 6-quarks these terms can contribute at the 0(5%) level for typical values of yc 
used in jet clustering. Therefore, the > 3-jet rate in heavy-quark events is expected 
to be modified relative to that in light-quark events both by the diminished phase- 
space for gluon emission due to the quark mass, as well as by kinematic effects in 
the definition of the jet clustering schemes. Such mass effects for jet rates have very 
recently been calculated |T^, |T3| at NLO in perturbative QCD|, and the quark-mass 
dependence can be expressed in terms of the running mass m}y{Mzo). The Aachen 
group has evaluated |^ the terms A^, and for massive 6-quarks at our preferred 
values of yc] these are hsted in Table 3. 

For illustration, the measured ratios R^/ R^'^'^ and R\/R^'^^, are shown in Fig. ^(a). 
R\/ R^'^'^ lies above unity for the E, EO, P and PO algorithms, and below unity for the 
D and G algorithms; note that all six data points are highly correlated with each other, 
so that the differences between algorithms are more significant than naively implied by 
the statistical errors displayed. For comparison, the corresponding QCD calculations 
of i^g/i^g'^'^ are also shown in Fig. ^(a), under the assumption of a flavor-independent 
strong coupling with an input value of as{Mz) = 0.118, for mb(M^o) = 3.0 ± 0.5 
GeV/c^. Under this assumption the calculations are in good agreement with the data, 
and the data clearly demonstrate the effects of the non-zero 6-quark mass, which are 
larger than the statistical error. For the translation from R\/ R^'^'^ to a^./a^'^'^ we used 
a value of the running 6-quark mass mb{Mzo) = 3.0 GeV/c^. 

For c-quarks mass effects are expected to be 0(1%) or less |^, which is much 
smaller than our statistical error of roughly 4% on R'^/R^'^'^ . The effects of non-zero 
c-quark mass, and of the light-quark masses, will hence be neglected here. We used 
values of A""^', B""^' and C"^' from Ref. [EBi. 



§In our previous study only the relevant tree-level calculations for 3-jet and 4-jet final-states 
were available. 
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Eqns. (2) were solved to obtain the ratios a^/a^'^* and a^/ a^'^^ for each jet algorithm. 
These ratios are listed in Table 2, together with the corresponding statistical and 
experimental systematic errors, and the hadronisation uncertainties. We then evaluated 
sources of uncertainty in this translation procedure. From an operational point of view 
these affect the values of the coefficients A, B and C used for the translation. For 
each variation considered the relevant A, B oi C were reevaluated, the ratios a\/a^'^'^ 
were rederived, and the deviation w.r.t. the central value was assigned as a systematic 
uncertainty. 

We considered a variation of ±0.5 GeV/c^ about the central value of the running 
6-quark mass mf,(M^o) = 3.0GeV/c^. This corresponds to the range 3.62 < mhlmh) < 
5.06 GeV/c^ and covers generously the values determined from the T system 



using QCD sum rules, 4.13 ± 0.06 GeV/c^, as well as using lattice QCD, 4.15 ± 0.20 
GeV/c^. It is also consistent with the recent DELPHI measurement of the running 
mass: mfe(M^o) = 2.67 ± 0.25{stat.) ± 0.34(/ra^.) ± 0.27(theo.) GeV/c^ g. The 
numerical accuracy on the coefficients A, B, and C is in all cases negligibly small on 
the scale of the experimental statistical errors. 

We considered the effects of the uncalculated higher-order terms in Eq. (2). In these 
ratios the effects of such higher-order contributions will tend to cancel. Nevertheless we 
have attempted to evaluate the residual uncertainty due to these contributions. We first 
considered 3-jet contributions and varied the NLO coefficient B; for each jet algorithm 
we varied simultaneously the renormalisation scale fi and a'^'^'^ in the ranges allowed 
by fits to the flavor-inclusive differential 2-jet rate ||15|[1 In addition, we considered 
next-to-leading order (NLO) contributions to the 4-jet rate. Although these enter 
formally at 0{a^) in Eq. (2), operationally they may be estimated by variation of the 
LO coefficient C*. Since the 4-jet rate has been calculated recently complete at NLO 
for massless quarks these terms can be estimated reliably. For our jet algorithms 
^ Heavy-quark mass and possible flavor-dependent effects are negligible on the scale of the large 
errors considered on a"*^* for this purpose. 
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and Uc values Dixon has evaluated the LO and NLO 4-jet contributions |3^. Based on 



these calculations we varied the coefficient C by ±100%. For each jet algorithm, at 
the chosen Uc value, the measured contribution to -R3 from >5-jet states was smaller 
than 1% and the corresponding 0{al) contributions to Eq. (2) were neglected. 

These uncertainties are summarised in Table The deviations for each varia- 
tion considered were added in quadrature to define a total translation uncertainty on 
a'^Jaf' and a\/af', listed in Table 2. 

7. Comparison of ag Ratios 

The Og/a"'^'' and a^/a^'^^ ratios are summarised in Fig. ^d. It can be seen that the 
ratios determined using the different jet algorithms are in good agreement with one 
another. 

For each jet algorithm ra, the statistical and experimental systematic errors were 
added in quadrature with the hadronisation and translation uncertainties (Table 2) 
to define a total error cr^ on = c or 6). For each flavor a single value of 

Og/a"'^'^ was then defined by taking the weighted average of the results over the six jet 
algorithms: 

«V«f^ = E<K/«:'ln (3) 

n 

where w\ is the weight for each algorithm: 

En l/< 

The average statistical and experimental systematic errors were each computed from: 



= Jj2^nm<^ln (5) 
y nm 

where E"^ is the 6x6 covariant matrix with elements: 

Km = ^n^m (6) 
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and 100% correlation was conservatively assumed among algorithms. The average 
translation and hadronisation uncertainties were calculated in a similar fashion. We 
then calculated the r.m.s. deviation on a^/a"''*' and a^/a"'^*, shown in Tabic 2, and 
assigned this scatter between the results from different algorithms as an additional 
theoretical uncertainty. The average translation and hadronisation uncertainties were 
added in quadrature together with the r.m.s. deviation to define the total theoretical 
uncertainty. 
We obtained: 

a^/a-'^ = 1.036 ±0.043(stat.)+°:°^^(syst.)l°:°^°(t/ieor|/) 
a\/af' = imA±QmS{stat)t'^^{syst.)t'^^{theory). 

The theoretical uncertainties are only slightly smaller than the respective experimental 
systematic errors, and comprise roughly equal contributions from the hadronisation 

and translation uncertainties, as well as from the r.m.s. deviation over the six jet 
algorithms. 

8. Cross-checks 

We performed a number of cross-checks on these results. First, we varied the event 
selection requirements. The thrust-axis containment cut was varied in the range 0.65 < 
I cos^tI < 0.75, the minimum number of charged tracks required was increased from 7 
to 8, and the total charged-track energy requirement was increased from 20 to 22 GeV. 
In each case results consistent with the standard selection were obtained. 

Next, wc included in the unfolding procedure (Eq. (1) and Section 4) the 'untagged' 
event sample, subsamplc 4 (Section 3), whose flavor composition is similar to the 
natural composition in flavor-inclusive Z° decay events, and repeated the analysis to 
derive new values of al/d^'^^ and a\/d^'^^. In addition, we repeated the unfolding and. 
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instead of fixing them to Standard Model values (Table 1), allowed the —>■ cc and 
—>■ bb branching fractions to float in the fit of Eq. (1). In both cases results 
consistent with the standard procedure were obtained [pO|] . 

We also considered variations of the flavor-tagging scheme based on reconstructed 
secondary vertices. In each case we repeated the analysis described in Sections 4-7 
and derived new values of a^a^'^'^ and a\/a^'^'^. Firstly, we used more efficient tags for 
primary bb and cc events. We applied the scheme described in Section 3, but with a 
looser definition of region (A) to include vertices with M^t^ > 1-8 or Pyt^ + ^O < ^^M^t^. 
We also removed the cut on the vertex impact parameter, S^tx, used to define the c- 
tagged sample, and region (B) was redefined to comprise only events with N^ig > 1 
and containing a SV with Pyt^ > 5 © P^tx + 10 > IbMyt^. Second, we repeated this 
modified scheme, but increased the efficiency for light-quark tagging by requiring tracks 
that miss the IP by at least Scxd to be counted in Ngig for the definition of the wrfs-tagged 
sample. Third, we did not use vertex momentum information for the tag definitions; 
we used instead only vertex mass information to define region (A): M^t^ > 1.8, and 
Region (B): M^t^ < 1.8, with the wrfs-tagged sample defined as in Section 4. Finally, 
we tried a variation in which we used event hemispheres as a basis for flavor-tagging, 
rather than jets as deflned in Section 3; this tag is similar to that used in our recent 
study of the branching fraction for Z^ bb [^. In all cases results statistically 
consistent with our standard analysis were obtained PD . 



We also performed an analysis using a similar flavor-tagging technique to that 



reported in our previous publication [11]. We counted the number of tracks per event, 
Ngig, that miss the IP by d > 30"^^. This distribution is shown in Fig. 7; the data are 
well described by our Monte Carlo simulation. For the simulation, the contributions 
of events of different quark flavors are shown separately. The leftmost bin contains 
predominantly events containing primary u, d, or s quarks, while the rightmost bins 
contain a pure sample of events containing primary b quarks, the event sample was 
divided accordingly into flve subsamples according to the number of 'signiflcant' tracks: 
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(i) Nsig = 0, (ii) Ns^g = 1, (iii) Ns^g = 2, (iv) Nsig = 3, and (v) A^..^ > 4. A similar 
formalism to that defined by Eq. (1) was applied using 5x3 matrices e and yielded 
values of Rf'/Rf, Rl/ Rf and R^/Rf consistent with those obtained in Sections 4 
and 5, but with larger statistical and systematic errors. Furthermore, we also applied a 
simpler version of this technique in which subsamples (ii), (iii) and (iv) were combined 
into a single c-tagged sample and a 3 x 3 flavor unfolding was performed. Again, this 
yielded values of R^'^^ / Rf\ R^/ Rf- and R\/ Rf^ consistent with those obtained in 
Sections 4 and 5, but with larger statistical and systematic errors p]|]. 

9. Summary and Discussion 

We have used hadron lifetime and mass information to separate hadronic decays 
into tagged 66, cc and light-quark event samples with high efficiency and purity, and 
small bias against events containing hard gluon radiation. From a comparison of the 
rates of multijet events in these samples, we obtained: 

al/af' = 1.036 ±0.043(stat.)i°:°^^(s?/st.)tn?^(t/ieory) 
al/af' = 1.004 ±0.018(stot.)l°;°^6^^^g^_^+o.oi8(^^g^^^^)_ 

We find that the strong coupling is independent of quark flavor within our sensitivity. 

For comparison with our previous result and with other experiments one can dis- 
cuss the test of flavor-independence in terms of the ratios a^'^'^ /a.f\ a^^/af'' and a^a"", 
although these quantities, by construction, are not independent of each another. We 
performed a similar analysis to that described in Sections 6 and 7 using, instead of 
Rl/Rf' and Rl/Rf', our measured values of Rf'/Rf, Rl/Rf and R\/Rf (Sec- 
tion 4) as a starting point. We obtained: 

«f7af = 0.987 ±0.010(stot.)l|]:[Jj2(s^,t.)+0-009(t/,eory) 
al/af = 1.023 ±0.034(stat.)in36(^2/^^-)-o:ou(i/^eor|/) 
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allaf = 0.993 ±0.016(stot.)l[i:lii(s2/st.)tn^?(t/ieor|/). 



These results are consistent with, and supersede, our previous measurements |[TT[], and 
are substantially more precise; they are also consistent with measurements performed at 
LEP using different flavor-tagging techniques H, ||, |lOl- A summary of these results 
is given in Fig. 8. Our comprehensive study, involving six jet-finding algorithms, and 
the inclusion of the resulting r.m.s. deviations of results as additional uncertainties, 
represents a conservative procedure. 
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Table 1: Compilation of the systematic errors for the EO algorithm and i/cut = 0.02. 
The first column shows the error source, the second column the central value used, 
and the third column the variation considered. The remaining columns show the 
corresponding errors on the values of R^'^'^ and R\/R^'''^] '+' ('— ') denotes the 
error corresponding to the relevant positive (negative) parameter variation. 



Source 


Center 


Variation 








Value 




+ 




-f 




tracking efSciencv 


correction 


off 


0.0020 




-0.0110 




2D imp. par. res. 


smear 


off 


-0.0100 




0.0080 




z track resolution 


smear 


off 


0.0010 




0.0120 




MC statistics 


0.8M 




0.0190 


-0.0190 


0.0091 


-0.0091 


B decav <^ 


5.51 trks 


±0.35 trks 


-0 0030 


-0 0026 


0.0135 


-0.0132 




0.697 


±0.008 


<0 0001 


0.0004 


0.0172 


-0.0191 


B frcLgin. slicipc 


Peterson 


Bowler 


0.0021 




-0.0216 




B meson lifetime 


1.56ps 


±0.05 ps 


-0.0021 


0.0022 


-0.0011 


0009 


B baryon lifetime 


l.lOps 


±0.08 ps 


-0 0003 


0003 


<0 0001 


-0 0000 


B baryoii prod. 


7.6% 


±3.2% 


0.0014 


-0.0016 


0.0021 


-0.0023 


B + X fraction 


0.192 


±0.05 


0.0011 


-0.0012 


-0.0013 


-0.0008 




0.2156 


±0.0017 


0.0022 


-0.0021 


0.0014 


-0.0014 




0.172 


±0.010 


0.0272 


-0.0294 


0.0044 


-0.0042 


C~1 ■FrfiP'TTi <' T ^ 


0.483 


±0.008 


0.0213 


-0.0211 


0002 


-0 0002 




T^etpTSOTi 


Rowler 


0.0042 




0006 






2 54 trks 


±0 06 trks 


0.0044 


-0.0048 


0006 


-0 0006 


D Herav <' t?,^*, 


2 48 trks 


±0 06 trks 


n 006Q 


-0.0074 


0.0012 


-0.0013 


decav <" Ti.^u '> 


2.62 trks 


±0.31 trks 


0039 


-0.0040 


-0.0004 


0003 


D° lifetime 


418 OS 


±0 004 ns 


-0.0001 


0.0001 


-0 0002 


0.0001 


D""" liflptimp 


1 054 ns 


±0 015 ns 


0.0001 


-0.0001 


-0.0001 


0.0001 


lifetime 


0.466 ps 


±0.017 ps 


0.0001 


-0.0001 


-0.0003 


0.0003 


Z)0 ^ ifO mult. 


0.402 


±0.059 


0.0088 


-0.0089 


0.0026 


-0.0026 


D+ K° mult. 


0.644 


±0.078 


0.0102 


-0.0120 


0.0027 


-0.0027 


Ds K° mult. 


0.382 


±0.057 


0.0012 


-0.0013 


0.0003 


-0.0003 


D° no TT*^ fraction 


0.370 


±0.037 


0.0069 


-0.0075 


0.0034 


-0.0034 


Z)+ no tt" fraction 


0.496 


±0.050 


0.0017 


-0.0018 


0.0029 


-0.0029 


Ds — !■ no tt'^ fraction 


0.348 


±0.035 


-0.0002 


0.0001 


-0.0003 


0.0003 


cc ^ D"*" + X fraction 


0.259 


±0.028 


0.0029 


-0.0034 


0.0001 


-0.0002 


cc ^ Ds + X fraction 


0.113 


±0.037 


-0.0025 


0.0019 


0.0002 


-0.0002 


cc ^ Ac + X fraction 


0.074 


±0.029 


-0.0051 


0.0044 


-0.0001 


-0.0001 


Ac decay < rich > 


2.79 


±0.45 trks 


0.0003 


-0.0002 


0.0024 


-0.0024 


Ac lifetime 


0.216ps 


±0.011 ps 


-0.0037 


0.0011 


-0.0006 


0.0001 


g — > 66 rate 


0.31 


±0.11% 


0.0001 


-0.0001 


-0.0038 


0.0039 


g ^ cc rate 


2.38 


±0.48% 


-0.0019 


0.0020 


-0.0015 


0.0016 


iiT" prodn. 


0.658trks 


±0.050 trks 


-0.0051 


0.0045 


-0.0061 


0.0058 


A prodn. 


0.124trks 


±0.008 trks 


-0.0007 


0.0009 


-0.0008 


0.0009 


Total Exp. Syst. 






0.0440 


-0.0480 


0.0300 


-0.0370 


Qo 


1 GeV 


ti.s GeV 


0.0074 


-0.0027 


0.0062 


-0.0237 




0.39 GeV 


±0.04 GeV 


0.0042 


-0.0008 


0.0015 


0.0012 


hadronisation model 


JETSET7.4 


HERWIG5.9 


0.0123 




-0.0383 




Total Hadronisation 






0.0150 


-0.0028 


0.0065 


-0.0450 
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Table 2: R\/ R^f-^ and a\/a^'^^ values and errors. 



Algorithm 

Vc 


E 
0.040 


EO 

0.020 


P 

0.020 


PO 

0.015 


D 
0.010 


G 
0.080 








central val. 


1.043 


1.066 


1.004 


1.058 


1.038 


1.040 




Stat. 


0.061 


0.016 


0.016 


0.010 


0.062 


0.086 




exp. syst. 


+0.065 
-0.075 


+0.044 
-0.048 


+0.046 
-0.046 


+0.039 
-0.039 


+0.062 
-0.067 


+0.074 
-0.085 




hadronisation 


+0.012 
-0.001 


+0.015 
-0.003 


+0.014 
-0.004 


+0.015 
-0.003 


+0.006 
-0.003 


+0.008 
-0.004 




total. 


+0.092 
-0.099 


+0.065 
-0.067 


+0.067 
-0.065 


+0.058 
-0.056 


+0.088 
-0.091 


+0.114 
-0.121 






L . ILL. O 


central val. 


1.031 


1.054 


1.004 


1.052 


1.032 


1.035 


01 7 


stat. 


0.046 


0.037 


0.041 


0.035 


0.051 


0.074 




exp. syst. 


+0.047 
-0.054 


+0.036 
-0.039 


+0.041 
-0.041 


+0.035 
-0.035 


+0.051 
-0.055 


+0.064 
-0.073 




hadronisation 


+0.009 
-0.001 


+0.012 
-0.002 


+0.012 
-0.001 


+0.013 
-0.003 


+0.005 
-0.002 


+0.007 
-0.001 




translation 


+0.001 
-0.002 


+0.005 
-0.006 


<0.001 


±0.008 


+0.003 
-0.005 


+0.004 
-0.006 




Rl/Rl<^' 




central val. 


1.050 


1.054 


1.048 


1.055 


0.964 


0.995 




stat. 


0.026 


0.019 


0.019 


0.017 


0.023 


0.032 




exp. syst. 


+0.038 
-0.042 


+0.030 
-0.037 


+0.027 
-0.037 


+0.028 
-0.035 


+0.038 
-0.041 


+0.035 
-0.036 




hadronisation 


+0.011 
-0.046 


+0.007 
-0.045 


+0.002 
-0.026 


+0.007 
-0.037 


+0.001 
-0.006 


+0.020 
-0.008 




total. 


+0.047 
-0.067 


+0.036 
-0.061 


+0.033 
-0.049 


+0.033 
-0.054 


+0.044 
-0.047 


+0.051 
-0.049 






r.m.s 


central val. 


0.989 


0.995 


1.018 


1.014 


1.009 


0.993 


0.011 


stat. 


0.018 


0.015 


0.017 


0.015 


0.021 


0.027 




exp. syst. 


+0.026 
-0.029 


+0.023 
-0.028 


+0.023 
-0.032 


+0.024 
-0.030 


+0.034 
-0.036 


+0.030 
-0.031 




hadronisation 


+0.008 
-0.032 


+0.005 
-0.034 


+0.002 
-0.022 


+0.006 
-0.032 


+0.001 
-0.005 


+0.017 
-0.007 




translation 


+0.016 
-0.015 


+0.013 
-0.015 


+0.011 
-0.014 


+0.017 
-0.018 


±0.012 


+0.008 
-0.009 
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Table 3: The coefficients A'^, B^, for tfie Next-to-Leading-Order calculation for mas- 
sive quarks. The numbers in parentheses represent the estimated numerical precision. 
Theoretical uncertainties in the computation of the coefficients derive from the 
'slicing parameter' used to isolate singular regions of phase space, as well as from the 
conversion to the MS quark mass parameter Effects from higher-order perturbative 
QCD contributions are discussed in the text. 







for mb{Mzo) (GcV/c^) = 


for mb{Mzo) (GcV/c^) = 


Algorithm 


Vc 


2.5 


3 


3.5 


2.5 


3 


3.5 


E 


0.040 


14.392(1) 


14.459(1) 


14.543(1) 


443(4) 


466(4) 


487(4) 


EO 


0.020 


24.850(2) 


25.024(2) 


25.231(2) 


277(4) 


291(4) 


310(4) 


P 


0.020 


24.850(2) 


25.024(2) 


25.231(2) 


63(4) 


67(4) 


75(4) 


PO 


0.015 


30.054(2) 


30.315(2) 


30.631(2) 


2(4) 


14(4) 


29(4) 


D 


0.010 


15.355(2) 


15.213(2) 


15.060(2) 


105(4) 


102(4) 


99(4) 


G 


0.080 


11.493(1) 


11.435(1) 


11.365(1) 


61(4) 


58(4) 


57(4) 







for mb{Mzo) (GcV/c^) = 


Algorithm 


Vc 


2.5 


3 


3.5 


E 


0.040 


27.91(1) 


28.27(1) 


28.71(1) 


EO 


0.020 


125.39(7) 


127.34(7) 


129.55(8) 


P 


0.020 


125.39(7) 


127.34(7) 


129.55(8) 


PO 


0.015 


202.8(1) 


206.1(1) 


209.4(1) 


D 


0.010 


84.30(6) 


82.83(6) 


81.19(6) 


G 


0.080 


65.55(4) 


64.60(3) 


63.56(3) 
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Table 4: Summary of translation uncertainties on the ag ratios for each algorithm; 
'+' (' — ') denotes the error corresponding to the relevant positive (negative) parameter 
variation. 



source 


Center 


Variation 


da's/a^"' (%) 
± 


fc^/ar (%) 

+ 




E-algo {yc = 0.04) 


mb{Mz) 


3.0GeV 


±0.5 


0.000 0.000 


-0.014 0.015 


H,as dep. 






-0.002 0.001 


0.006 -0.006 


^ 4jet contrib. 


C 


±C 


<0.001 <0.001 


-0.001 0.001 


lotai 






0.001 -0.002 


0.016 -0.015 




EO-algo {jjc = 0.02) 


mb{Mz) 


3.0GeV 


±0.5 


0.000 0.000 


-0.014 0.012 


H, tts dep. 






-0.005 0.004 


0.005 -0.005 


> 4jet contrib. 


C 


±C 


-0.002 0.003 


-0.001 0.002 


lotai 






0.005 -0.006 


0.013 -0.015 




P-algo (yc = 0.02) 


mb{Mz) 


3.0GeV 


±0.5 


0.000 0.000 


-0.012 0.009 


11, as dep. 






<0.001 <0.001 


-0.002 0.002 


> 4jet contrib. 


C 


±C 


<0.001 <0.001 


-0.005 0.007 


Total 






<0.001 <0.001 


0.011 -0.014 




PO-algo {yc = 0.015) 


mb{Mz) 


S.OGeV 


±0.5 


0.000 0.000 


-0.017 0.015 


^, as dep. 






-0.007 0.005 


-0.001 <0.000 


> 4jet contrib. 


C 


±C 


-0.004 0.006 


-0.006 0.008 


Total 






0.008 -0.008 


0.017 -0.018 




D-algo {yc = 0.010) 


mb{Mz) 


S.OGeV 


±0.5 


0.000 0.000 


0.011 -0.010 


H, as dep. 






-0.005 0.002 


-0.005 0.003 


> 4jet contrib. 


C 


±C 


-0.002 0.002 


-0.003 0.003 


Total 






0.003 -0.005 


0.012 -0.012 




G-algo {ye = 0.08) 


mb{Mz) 


S.OGeV 


±0.5 


0.000 0.000 


0.010 -0.009 


/i, as dcp. 






-0.005 0.003 


0.005 -0.003 


> 4jet contrib. 


C 


±C 


-0.002 0.003 


0.001 -0.001 


Total 






0.004 -0.006 


0.008 -0.009 
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Figure 1: The distributions of (a) the vertex mass, M^t^, and (b) the vertex momentum, 
Pvtx, in our data sample (points); the simulated distributions are shown as a histogram 
in which the contributions from events of different primary quark flavor are indicated. 
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Figure 2: The two-dimensional distribution of vertex momentum P^jtx vs. vertex mass 
Mytx (see text), (a) Data; (b) all-flavors simulation; (c) bb event simulation; (d) cc event 
simulation; (e) qiqi simulation. In (f) the regions used for 6-tagging (A), c-tagging (B) 
and no-tagging (C) are indicated (see text). 
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Figure 3: The distribution of vertex impact parameter, 5ytx, for events containing 
vertices in region (B): data (points); the simulated distribution is shown as a histogram 
in which the contributions from events of different primary quark flavor are indicated. 
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Figure 4: The distribution of the number of tracks per event that miss the IP by at 
least 2(7 in terms of their impact parameter in the plane normal to the beamline, in 
events that contain no reconstructed vertex (see text); data (points). The simulated 
distribution is shown as a histogram in which the contributions from events of different 
primary quark flavor are indicated. 



34 



1.6 
1.4 



1.2 



1 — I I I I I 1 1 — I — I — I I I I I 

Subsamplej: uds c b No-tag 
Data { f { A 
MC 



— ■ CO 

DC 



1.0 
0.8 
1.4 



= 1-2 

OJ CO 



- CO 1 .0 



— ^ - - .A. 3 

(a) Raw data 



J I I L_l_ 



J I I I I I L_l_ 



SLD 



Sample i: uds c b 

} \ \ 



0.8 



0.6 



(b) Unfolded 

J I I I I I 



J I I I I I LJ_ 



4-98 
8358A11 



10" 



10" 



Figure 5: (a) The raw measured ratios -R3/-R3", 1 < J < 4, vs. yc for the 4 subsamples 
(see text); data (points with error bars), and simulation (hues joining values at the 
same yc values as the data), (b) The unfolded ratios R\/Rf^, i — h,c,uds, vs. yc for 
the 3 primary event flavor groups. Only statistical errors are shown. In (b) points 
corresponding to a common y^. value have been displaced horizontally for clarity. 
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Stat, error only 




Figure 6: (a) The measured ratios R\/R^'^^, and (b) the corresponding translated 
ratios a\/a^'^^ {i — c,h). The arrows in (a) indicate the range of the theoretical 
prediction described in the text for values of the 6-quark mass in the range 2.5 < 
miy{Mza) < 3.5 GcV/c^, with the arrow pointing towards the lower mass value. In (b) 
The weighted average over the six algorithms is also shown. In all cases only statistical 
error bars are displayed. 
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Figure 7: The distribution of the number of tracks that miss the IP by at least 3a in 
terms of their impact parameter in the plane normal to the beamline (see text): data 
(points); the simulated distribution is shown as a histogram in which the contributions 
from events of different primary quark flavor are indicated. 
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Figure 8: Summary of measurements of al/af' {i = uds, c or b) from experiments at 
the resonance. We derived the ALEPH a'^'^^/af^ value from their measured value 
of a'^^^/a'^f, as well as the five bracketed LEP values of ct^/o;"" from the measured 
values of a^/a"'^*"^, by assuming a"" = J2uds,c,b P ^^l^ where /* is the Standard Model 
branching fraction for Z° decays to quark flavor i. 
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